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ABSTRACT 

Recent theoretical and observational studies both argue that the merging 
of double carbon-oxygen white dwarfs (WDs) is responsible for at least some 
Type la supernovae (SNe la). Previous (standard) studies of the anticipated SN 
birthrate from this channel have assumed that the merger process is conservative 
and that the primary criterion for explosion is that the merged mass exceeds the 
Chandrasekhar mass. Han & Webbink (1999) demonstrated that mass trans- 
fer and merger in close double WDs will in many cases be non-conservative. 
Pakmor et al. (2011) further suggested that the merger process should be vi- 
olent in order to initiate an explosion. We have therefore investigated how the 
SN la birthrate from the double-degenerate (DD) channel is affected by these con- 
straints. Using the binary-star population-synthesis method, we have calculated 
the DD SN la birthrate under conservative and non-conservative approximations, 
and including lower mass and mass-ratio limits indicated by recent smoothed- 
particle-hydrodynamic calculations. The predicted DD SN la rate is significantly 
reduced by all of these constraints. With dynamical mass loss alone (violent 
merger) the birthrate is reduced to 56% of the conservative rate. Requiring that 
the mass ratio q > 2/3 further reduces the birthrate to 18% that of the standard 
assumption. An upper limit of 0.0061 SNuM, or a Galactic rate of 4.6 x 10 yr -1 , 
might be realistic. 

Subject headings: binaries: close — stars: evolution — supernovae: general 
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Introduction 



The double-degenerate (DD) merg er was suggested as a possibl e channel for type la su- 



pernovae ( SNe l a) in the early 1980's ( ITutukov fc Yungelson Ill98ll ; llben fc Tutukov 111984 



Webbink 1 1 19841 ) . In this DD model, two carbon-oxygen (CO) white dwarfs (WDs) can 



produce a SN la while merging if their total mass is larger than the Chandrasekhar mass 
(M c h). The DD model can naturally explain the lack of H and He emission in the spectra 
of SNe la and some super-luminous SNe la, but it has a major difficulty in explaining the 
similarities of most SNe la since the merger mass has a relatively wide range, ~ 1.4 — 2.OM 
rtWang et al. I lioiolfl Meanwhile, alth ough the SNe la rate predicted from the DD model 



is comparable to that of observations (lYungelson et al. 



200ll ; iRuiter. Belczvnski fc Frv"e7l booguWang et al. I hold : lYu fc Jeffervl bold ), theoreti- 



199 



Han 



19981; Nelemans et al. 



cal studies suggest that the merger of two WDs is more likely to lead to an accretion- 



induced collapse (AIC) to form a neutron star ( INomoto &: Iben Ill985l ; ISaio &: Nomoto 111985 



Timmes. Woosley fc Taam 1 119941 ) . Consequently, the DD model has been disadvanta ged in 
comparison to the single-degenerate (SD) model for a long time (IWang fc Han 1 120121 ). 



Piersanti et al. I ( 120031 ) suggested that, under the right conditions, the DD merger 



process could be quite violent and lea d to a SN la explos ion rather than AIC (see also 
Yoon. Podsiadlowski fc Rosswog I ( 120071 ) ; IShen et al. I (120121) ). Fully three-di mensional sim- 



ulations of a violent merger of two CO WDs by iPakmor et al. I ( 120101 . 120111 ) show that the 
merger of two equal-mass CO WDs (~ 0.9M^) can explain the form ation of sub-luminous 



1991bg-like events (see also Ivan Kerkwijk. Chang fc Justham I (120101 ) ). The IPakmor et al. 



( 120121 ) simulation of a DD merger with masses of 1.1M and O.9M shows good agreement 
with properties of normal SNe la. Observat ionally, some known or potential DD systems suc h 



as WD 2020-425 (INapiwotzki et al. 



SBS1150+599A (ITovmassian et al. 



20071 ). V458 Vufoe culae (IRodriguez-Gil et al. 1 120101 ) 



2010h and GD687 JGeier et al. I boiol ) represent good 



candidates for SN la progenitors since they have total masses close to M c h and orbital pe- 
riods short enough to merge within a Hubble time. Furthermore, extensive searches have 
found no surviving compan ion to the T ype la supernovae respo nsible for SNR 0509-67.5 
( ISchaefer fc Pagnotta 1120121 ) or SN 1572 (IKerzendorf et al. 1120121 ). whence it is argued that 
the progenitors must have been DD systems. Hence, it appears that at least some, if not 
most, SNe la come from DD mergers. 



Previous studies for the DD model generally assume that a SN la is produced if the 



1 It can be argued, however, that in the DD scenario it is the primary (exploding) WD which matters 
most, not the tota l mass, since it is the deton ation in the exploding WD which is primarily responsible for 
synthesizing 56 Ni (jPakmor et al. Iboid Eoij) . 
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total mass of a DD system is larger than M c h. This assumption implies that the DD merger 
process is always conservative. However, if accretion is spherically symmetric, mass transfer 
can be dramatically non-conser vative both during stable mass transfer and in a violent 
merger. lHan fc Webbink I (119991 ) (HW) studied the stability and energetics of mass transfer 



in double WDs. They showed that the expelled mass fraction of an interacting double WD 
depends on the component masses and can be as large as 50% of the donor for a system with 
masses of lM ffl and Q.5M&. Meanwhile, smoothed-particle-hy drodynamic (SPH) simulations 
(jYoon. Podsiadlowski fc Rosswog 1120071 ; iPakmor et al. Il201ll ) have shown th at not all double 



CO W Ds with a total mass larger than M c h can produce SNe la. For example, IPakmor et al 



(120 111 ) find that only double CO WDs in which the more massive white dwarf exceeds 
~ O.9M and the mass ratio roughly exceeds 0.8 robustly reach the conditions required for 
initiating a detonation, leading to a SN la explosion. Obviously, some additional constraints 
should be placed on the DD model for triggering SNe Icfl. 

In sect. 2 we investigate the merger mass of double WDs including non-conservative 
mass transfer, and we examine the impact of these constraints and those arising from SPH 
simulations on the SNe la birthrate in sect. 3. Conclusions are drawn in sect. 4. 



2. The Merger Mass 

As the orbit of a double WD system decays due to the emission of gravitational radi- 
ation/waves, the less massive component will be the first to fill its Roche lobe and become 
the mass donor; the more massive component becomes the accretor. Unlike non-degenerate 
donor stars in other compact binaries, the WD donors in double WDs lie deeper within 
the potential well of the accreting star, resulting in important consequences for mass trans- 
fer. Assuming energy conservation and hence no energy sources other than gravitational 
potential energy are at play, HW analyzed the initial stability of this mass transfer and the 
fraction of the mass-transfer stream to be expelled from the binary. From this point of view, 
the mass transfer is conservative if the accretion luminosity is less than the Eddington limit. 
Any accretion luminosity in excess of the Eddington limit is assumed to be absorbed in the 
accretion flow. In terms of energetics, radiative losses represent accretion energy not used 
to power a mass outflow (HW: Eq (21)). The fraction of matter accreted by the primary is 
used to power the outflow and can be determined from HW: Eq (22). For component masses 



2 From the point of view of preventing an off-center C ignition, lYoon. Podsiadlowski fc Rosswog I (|2007l ) 
favored less massiv e double CO WDs to produce SNe la. Since an off-center C ignition is almost inevitable 



(|Shen et al. 1120121 ). and much of the p hysics for the therma l evolution of merger remnants is unclear, we 



only consider here the constraints from lPakmor et al. I (|201l[ ). 
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M\ (accretor) and M 2 (donor), we adopt HW: Eqs (17), (22) and (23) to calculate the mass 
transfer rate, M 2 , and the fraction of matter to be accreted by the primarj{^, (5. 

There are two regimes for computing the final total mass of the merger, M. If the 
merger process is dynamical, (3 can be applied as a single step process and M = Mi + (3M 2 . 
If mass transfer is stable, M is obtained by integrating over the mass-transfer epochQ, using 
M = M l + J f3(t)M 2 (t)dt. 



3. Type la supernovae from the DD model 



3.1. The formation of double CO WD systems 

To investigate the SN la birthrate due to DD mergers, we require a population of 
close double CO WD binarie s, to be formed from interactions in a population of primordial 
binary stars (e.g. lHan I (119981 ) ). We first perform a Monte-Carlo simula tion to obtain a stellar 



popu l ation , then evolve them in a rapid binary evolution code (RBEC) (IHurley. Tout fc Pols 
2000l . 120021 ) to obtain a sample of close CO WD systems. 



In the Monte Carlo simulation, all stars are assumed to be members of binaries and 



have circular orbits. The primaries follow the in itial mass function of iMiller &: Scalo I (ll979[ ) 
and are generated according to the formula of lEggleton. Tout fc Fitechett I ( 119891 ). in the 
mass range 0.08 to 100M Q . The secondary mass, also with a lower limit of O.O8M , is then 
obtained from a constant mass-ratio distribution. The distribution of orbital separations a 
is taken to be constant in loga for wide binaries. This separation distribution has been used 
in many Monte- Carlo simulations and implies an equal number of wide binary systems per 
logarithmic inter val and app roximately 50 per cent of stellar systems with orbital periods 
less than 100 yr ( Han Ill998l ). Long-orbital-period binaries are effectively single stars. 



The formation of close double CO WD systems depends significant ly on the cr itical mass 
ratio for dynamical instability q c and common envelope (CE) evolution ( IHan ll998h . A binary 
experiences stable Roche lobe overflow (RLOF) when q < q c and C E evolution otherwise, 
where q = donor/accretor. In this study we adopt equation (57) of IHurley. Tout &: Pols 



( 120021 ) for q c when the mass donor is on the first or asymptotic giant branch (AGB), and 



3 Note that the model of HW takes the accretion to be spherical. If mass transfer is through a disk, then 
excess energy can be radiated away from the poles and (3 may not be appropriate. 

The fitting formulae and a table for the merger mass from the two regimes for various component masses 
may be obtained by email from xuefeichen717@hotmail.com. 
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let q c = 4 when the mass donor is on the ma in sequence or in the Hertzsprung gap, as 



supported by detailed binary evolution studies ( IHan. Tout &: Eggleton 



2002 



2000 



Chen & Han 



200 3|). For CE evolution, we us e the standard energy formalism ( Ivan den Heuvel 



1976 



Webbink 111984 iLivio fe Soker Ill988l ). that is, the CE is ejected if a ce AE orh > E hin d, where 



Ai^orb is the orbital energy released, a ce is the CE ejection efficiency, and -Ebind is the binding 
energy of the envelope and can be written as G MdM env / (\Rd) (where M env and Rd are 
the mass, envelope mass and radius of the donor, respectively, and A is a structure parameter 
that depends on the evolutionary stage of the mass donor). We combine a ce and A into one 
free parameter, setting q re A = 1.5 to reprod uce the number of the DD objects in the Galaxy 
as in previous studies (IWang et al. Il2009bl ). 

RBEC distinguishes three types of WDs, namely helium (expected only in binaries), CO 
and oxygen- neon (ONe), respectively. The critical condition for CO and ONe WDs is the 
core mass at the base of AGB, -Mc,bagb- CO WDs are produced when M^bagb 

< 1.6M , 

and ONe WDs w hen 1.6M < M CiB agb < 2.5M Q , where M CjB agb is given in Eq (66) of 



Pols et al. I (119981 ). The mass of the CO core, and the final mass of CO WDs is determined 
by the L — M c relation, i.e. Eqs(37) and (39) of that paper. 



3.2. Constraints on the progenitors of SNe la 



There are two basic constraints on double CO WDs as the progenitors of SNe la in the 
conservative case: (i) Mi + M2 > M c Jf|, and (ii) the CO WDs are close enough to merge in 
a Hubble time. Th e timescale for two compo nents to merge by gravitational wave radiation, 
£gw> is written as (ILandau fc Lifshitz Ill97ll ) 



t, 



GW 



(M 1+ M^ 



(1) 



where P is the orbital period in hours, tew i n years and Mi, M2 in solar mass. 



In the non-conservative case, (i) should be redefined as the final mass M > M c h (see 
sect. 2). Furthermore, AIC should be avoided in order to success fully trigger a SN la . 



We therefor e need additional constraints from dynamical simulations. iPakmor et al. I ( 12010 



20 111 . |2012| ) showed that an initial detonation at the onset of C ignition can lead the DD 
merger to avoid AIC and explode, and that DDs with Mi > O.9M and a mass ratio 
M1/M2 > 0.8 can achieve conditions for initiating such a detonation. Their studies also 



Han & Podsiadlowski 


( 


2004): 


Mens: et al. ( 


2009); 


Wang et al. 



d2009ah . 
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showed that mergers with M\ ~ O.9M are very promising candidates for explaining sub- 
luminous SNe la, while a normal SN la can be produced from the double CO WDs with 
Mi = I.IMq. Since the SN la luminosity is determined by the amount of 56 Ni synthesized 
in the explosion, these simulation results seem to indicate that the primary mass is related 
to the final production of 56 Ni, and therefore determines whether a normal or sub-luminous 
SN la is produced. Meanwhile, we do see a lot of variation in SN la lightcurves, and thus 
one might expect to see a variety among exploding WD masses provided by the DD models. 
In the absence of other constraints, we study the SN la birthrate assuming two lower mass 
limits for the DD accretor, i.e. Mi > 0.9 or 1.0M Q , respectively. 

Another constraint from the SPH simulations is the mass ratio of the CO DDs, q = 
M 2 /Mi. It has a critical value, ~ 0.8, below which the merger is not violent enough to ignite 
a detonation (Pakmor et al. 2011)0. The relation between the violence of a merger and the 
mass ratio can be simply understood from Fig. 2 of HW. In the dynamical instability region, 
Cad — Cl < 0, where Cad and Cl are the adiabatic mass-radius exponent of the donor and the 
mass-radius exponent of its Roche radius, respectively. With the mass ratio closer to unity, 
the value of | Cad — Cl I becomes larger and larger. Since it is a negative number, the process 
of self- amplifying M 2 (Eq (1) in that paper) is thus faster and faster, resulting in a more and 
more violent merger process. We arbitrarily relax the critical mass ratio from 0.8 to 2/3 in 
our study to give an upper limit for the SN la birthrate from the SPH simulations. 

In brief, based on the sample of double CO WDs which can merge in a Hubble time, we 
compute the SN la birthrate for the following cases : 

(i) M = M 1 + (3M 2 > M ch ; 

(ii) M = M 1 + J (3(t)M 2 (t)dt > M ch ; 

(iii) As case (ii) with Mi > 0.9M Q . We choose case (ii) as the basic constraint since it 
includes a larger parameter space than case(i) (see Fig.l). The overall SNe la rate from case 
(ii) can then be considered as an upper limit from the non-conservative assumption. 

(iv) similar to (iii) but Mi > l.OM , and 

(v) an additional constraint, q = M 2 jM\ > 2/3, is placed on case (iii). 

Figure 1 shows these constraints in the Mi — M 2 plane (left panel), as well as CO DDs 
with a total mass Mi + M 2 > M c h which can merge in a Hubble time (right panel). To obtain 
these DDs, we generate a population of 10 7 binaries (with a total mass of ~ 1.26 x 10 M ) 



6 Pakmor et al. (2011) found that the critical value of q probably changes with primary WD mass; i.e. 
more massive primaries can merge with lower mass ratios. 
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as described in sect. 3.1 and evolve them in RBEC These DDs cover a wide range in the 
Mi — M 2 plane, and also a large range of j3. For q = 2/3, the /3 decreases from 0.7 to 0.54 
as Mi increases from 0.827 to 1.2M . The overall parameter spaces of SN la progenitors in 
cases (i)-(iv), i.e. to the right of the red lines, are reduced by different degrees relative to 
the conservative assumption, indicating corresponding decreases of the DD SN la rate under 
the non-conservative assumption and with dynamical constraints. 



3.3. Birth rates of SNe la 

The CO WDs above are used to study the SN la birthrates in various cases, which are 
finally normalized to (a) a single starburst of 1O 11 M and (b) a constant star formation 
rate of 5M Q yr _1 over the past 15 Gyr, to resemble our Galaxy. The results are shown in 
Fig.2. We see that, in the conservative case, the predicted SN la birthrate is 0.033 SNuM@ 
(corresponding to a Galactic birth rate of 2.5 x 10 _3 yr _1 ) for the constant star-formation 
at 15 Gyr (left panel). In the non-conservative case, the predicted overall rate is reduced 
to 0.019 SNuM (1.4 x KTV -1 ) and 0.024 SNuM (1.8 x lO^yr" 1 ) for cases (i) and (ii), 
respectively. If we introduce additional restrictions from dynamical simulations, the overall 
SN la rate is further reduced to 0.017 SNuM (1.3 x 10- 3 yr~ X ), 0.0099 SNuM (7.4 x 10- 4 yr -1 ) 
and 0.0061 SNuM (4.6 x 10~ 4 y r_1 ) f° r cases (iii)-(v), respectively. The birthrate in case (v) 
is only about 18% of that of the conservative case. Note that the SN la birthrate from case 
(v) might / could represent an upper limit based on constraints from recent SPH simulations 
(see sect. 3. 2). 

Cases (i)-(iv) present similar delay-time distributions to that of the conservative case, 
that is, a peak around 0.8 — 3.2 x 10 8 yr with a tail decaying as t~ x , where t is the delay time 
(Fig. 2, right). The shape of this distribution is consistent with both theoretical analyses and 
observations (see Maoz et al. 2011; Wang & Han 2012 and references therein). However, in 
case (v), one sees a small peak around logt(yr) = 8.7 and a 'dip' just before logt(yr) = 8.5. 
RBEC produces CO WDs nearly uniformly distributed in orbital period. The CE+CE 
channel produces extremely short orbital periods and the RLOF+CE channel gives relatively 
longer orbital periods, with some overlap between. Requiring q > 2/3 removes I0W-M2 
binaries, mainly from the RLOF+CE channel. These I0W-M2 binaries are produced in two 
ways; i) the AGB progenitors of M 2 before the CE phase have low masses and generally lower 
binding energies and produce double CO WDs with longer orbital periods for a given M%, ii) 
binaries with short orbital periods before CE ejection produce double CO WDs with short 



7 1 SNuM= lSN(100yr)- 1 (10 10 M Q )- 1 



0.2 0.4 0.6 0.8 1 1.2 1.4 ° 0.8 1 1.2 

If, «f„) U, (Jf ) 

Fig. 1. — Left panel: constraints on the progenitors of SNe la of double CO WDs in the 
Mi — M 2 plane. The black lines are contours of the ejected mass fraction 1 — with values 
of, from left to right, 0.0, 0.1, ... and 0.9, respectively. The red lines mark contours of merger 
mass Mi + M2 = M c h (conservative case, solid line), Mi + (3M2 = M c h (non-conservative, 
dotted line) and M l + j f3(t)M 2 (t)dt = M ch (non-conservative, dashed line). The vertical 
lines give the boundaries for Mi = 0.9 and l.OM . The areas to the right of each colour 
line represent M > M c h and hence DD SN la progenitor candidates under the corresponding 
constraints. Right panel: CO DDs from our population synthesis with Mi + M2 > M c h 
which can merge within a Hubble time, the blue dots have mass ratios q = M2/M1 > 2/3. 
The black dash-dot-dot-dotted line shows the boundary of q = 2/3, and the red lines are the 
same as those in the left-hand panel. 
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orbital periods. Hence the q > 2/3 constraint removes both the longest- and shortest-period 
binaries from the RLOF+CE channel, leaving a relatively larger number of progenitors with 
orbital periods of ~ 3 hrs upon emerging from the CE+CE channel (corresponding delay 
times of logt = 8.7). The 'dip' arises between the surviving CE+CE and RLOF+CE samples 
since short-orbital-period systems from the RLOF+CE channel have been removed. 



Conclusion 



We have shown that the DD model for SNe la is less effective than previously assumed 
if we consider the non-conservative nature of the DD merger process and constraints on 
triggering SN la explosions provided by recent SPH simulations. The overall SN la rate 
is reduced to a value rs 0.56 — 0.72 that of previous studies (or 0.018 — 0.024 SNuM) by 
introducing the non-conservative approximation, while it is likely below 0.18 of that value 
(or < 0.006 SNuM) if we also consider the SPH constraints. 

Theoretical estima tes of the SN la rate from the SD model show a wide range of values 
from w 0.001 SNuM jRuiter. Belczvnski fc Frver I l2009h to > 0.01 SNuM ( Han & Pod- 
siadlowski 2004) (see also Mennekens et al. 2010), which strongly depends o n the model 



assumptions and input parameters used in various studies. IWang et al. I ( 120101 ) performed a 



study using the same synthesis model and Galactic (constant) star formation rate as used in 
this study, and found a rate of 0.029 SNuM (2.15 x 10~ 3 yr _1 )- By direct comparison, the DD 
model appears unlikely to be better than the SD model in explaining the SN la birthrate. 

It was found by Pakmor et al. (2012) that an extended envelope enshrouding the merging 
WDs is not produced in violent merging events. However, we consider the possibility that a 
hot extended envelope may arise following a non-violent non-conservative DD merger. Since 
the disrupted material of the secondary must also contain most of the original orbital angular 
momentum, the outermost layers of the merger probably form a centrifugally supported disc. 
Our study therefore suggests that the lost or unaccreted material from the DD merger would 
take the form of a hot envelope plus disk , similar to that indicated by the SPH simulations of 
Yoon. Podsiadlowski fc Rosswog I ( 120071 ) in which the scale of the hot envelo pe is small (i. e. 



< 10 cm). The absence of early ultraviolet-optical emission in SN 201 lfe (INugent et al. 



201 ll ) is compatible with the DD model with a small extended envelope, while a disk would 
cause some continu um polarization such as that detected at red wavelengths in SN 201 lfe 
( ISmith et al. 1 1201 lh . although other causes of polarization cannot be ruled out. 
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Fig. 2. — Left panel: the evolution of birthrate of SNe la for a constant Pop I star formation 
rate (5M yr _1 ) from the DD model. Right panel: as the left panel for a single starburst of 



10 n M ( 



©■ 



